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ABSTRACT 

Performance  improvement  of  a  carbon  dioxide  gas-dynamic 
laser,  due  to  regeneration,  was  studied  theoretically.   The 
study  determined  to  what  extent  preheating  the  nitrogen, 
before  mixing  it  with  the  carbon  dixoide  and  water,  would 
effect  available  laser  output  energy.   Combustion  chamber 
pressure  and  nitrogen  inlet  temperature  were  incremented 
by  100  psia  and  3  00  K  respectively.   Due  to  the  fact  that 
little  dissociation  occurred  at  temperatures  studied,  no 
pressure  dependence  was  found  for  the  throat  or  combustion 
chamber  temperatures.   A  method  of  preheating  or  regenerating 
the  nitrogen  is  to  use  a  gas-to-gas  heat  exchanger  at  the 
exhaust  of  the  diffuser.   Small  signal  gain  was  improved 
from  0.3  5  to  0.7  8  m~  .   Limiting  power  extraction  was 
increased  by  42%.   Optimum  nitrogen  temperatures  were 
determined  to  be  13  50  K  for  small  signal  gain  and  2400  K 
for  optimum  power  extraction.   For  typical  operating  con- 
ditions the  fractional  pressure  loss  across  the  heat 
exchanger  was  found  to  be  0.1  atmospheres. 
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LIST OF  SYMBOLS 

2 

A     -   exchanger  total  heat  transfer  area  on  one  side.   (ft  ) 

2 

A     -   exchanger  free  flow  area.   (ft*) 

23     -1 
A     -   Avagadro's  Number  (6.02  x  10   Mole   ) 

C     -  constant  pressure  specific  heat 

f  -      Fanning  friction  factor 

2 

G     -   exchanger  flow  stream  mass  velocity  [lbs/(Hr-Ft  )] 

g     -  gravitational  constant  [32.2  (ft/sec  )  (jr— •)  ] 

-34  m 

h     -   Planck's  constant  [6.62  x  10    J-Sec] 

2 

I     -   saturation  intensity  (Kw/m  ) 

2 

I     -   output  beam  intensity  (Kw/m  ) 

J     -   rotational  quantum  number. 

-23 

k     -   Boltzrnann '  s  constant  [1.38  x  10    J/K] 

K     -   contraction  loss  coefficient  for  flow  at  heat 
exchanger  entrance. 

K     -  expansion  loss  coefficient  for  flow  at  heat 
exchanger  exit. 

M     -   Mach  number 

m     -  mass  of  one  molecule  of  xx. 

XX 

m     -  mass  flow  rate  of  xx. 

XX 

P*    -  pressure  in  cavity  (region  4  in  Figure  1)  psia 

q     -   rate  of  heat  transfer  from  hot  gases  to  cold 
nitrogen. 

R     -  universal  gas  constant 

TN    -   nitrogen  temperature  (K) 

T,     -   combustion  chamber  temperature  (K) 


T2    -  mixing  chamber  temperature  T2  =  T,  (K) 

T3    -   nozzle  throat  temperature  (K) 

T4    -  translational-rotational  temperature  in  the  cavity  (K) 

v,    -  upstream  specific  volume  of  gases  with  respect  to 
heat  exchanger. 

v2    -  downstream  specific  volume  of  gases  with  respect  to 
heat  exchanger. 


m 


2 


Vl  +  V2 


mean  specific  volume  =  _- 


a    -   small  signal  gain  coefficient  (m  ) 

y    -  C  /C   =1.4  for  this  mixture  of  gases 
p  v 

£    -  effective  optical  cross  section 
o 

a  -     ratio  of  free  flow  area  to  frontal  area 

c    -  optical  broadening  collisional  cross-section  of 
xx  (cm  ) 

8.    -   characteristic  vibration  temperature  of  CO   for 
1      v][.   (1995  K)  z 

0_    -   characteristic  vibration  temperature  of  C09  for 


V2.   (960  K) 


0-  -  characteristic  vibration  temperature  of  N2  (33  80  K) 

0  -  characteristic  rotational  temperature  of  C09 

r  (2.88  K) 

ty  -  mole  fraction  of  xx. 

p  -  density  of  gases 

13 

v  -  laser  frequency  (2,83  X  10   )  Hz 

v  -  optical  collision  frequency 

t  -  radiative  lifetime  (5.38  sec) 


u    -   reduced  mass  of  xx 
xx 


-6 


wavelength  of  laser  radiation  (10.6  X  10   m) 
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I.   INTRODUCTION 

The  carbon  dioxide  gas  dynamic  laser  has  been  considered 
a  prime  candidate  for  many  high  energy  laser  programs. 
One  problem  that  has  always  plagued  the  high  energy  laser 
program  is  that  of  excessive  weight  and  volume.   These 
presently  are  the  factors  that  are  impeding  the  development 
of  a  tactical  Airborne  laser  weapon  system.   A  heat  exchanger 
can  be  placed  at  the  exhaust  of  the  diffuser  to  preheat  the 
nitrogen  gas.   This  system  significantly  reduces  the  weight 
and  volume  of  a  given  gas  dynamic  system  while  maintaining 
a  constant  level  of  energy  output. 

The  combustion  chamber  properties  and  nozzle  throat 
conditions  were  determined  by  using  a  computer  program 
developed  by  Sanford  Gordon  and  Bonnie  McBride  at  the  NASA 
Lewis  Research  Center  in  Cleveland,  Ohio.   The  Program  and 
a  detailed  description  of  its  use  are  contained  in  reference 
[1] .   Although  the  program  was  written  with  a  rocket  engine 
analysis  in  mind,  it  is  adequate  for  computing  equilibrium 
and  frozen  composition  conditions  in  combustion  gases  used 
in  gas  dynamic  lasers. 

After  the  computer  determined  conditions  up  to  and 
including  the  nozzle  throat,  hand  calculations  were  made 
using  methods  specified  in  reference  [2] .   According  to 
calculations,  the  maximum  available  energy  can  be  increased 


by  a  factor  of  5  to  10.   The  small  signal  gain  as  a  function 
of  nitrogen  temperature  reached  a  maximum  value  of  0.773. 
This  maximum  indicates  a  best  operating  temperature. 

The  diffuser  is  used  to  recover  the  pressure  in  the  flow 
stream  and  must  be  matched  to  the  heat  exchanger. 

A  gas  to  gas  heat  exchanger  is  used  to  insure  that  the 
proper  temperature  difference  is  attained  from  the  cold 
nitrogen  side  to  the  hot  nitrogen  side  (see  Figure  1) .   The 
pressure  must  not  be  allowed  to  decrease  to  atmospheric 
pressure  at  mean  sea  level.   This  can  be  done  by  use  of 
procedures  outlined  in  references  [3],  [4],  and  [5]. 
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Figure    1         BASIC    ARRANGEMENT    FOR   A    REGENERATIVE    GAS 
DYNAMIC    LASER   WITH    MOPA   CONFIGURATION. 
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II.   THEORY 

A.   BASIC  THEORY  OF  A  GAS  DYNAMIC  LASER  AMPLIFIER 

Carbon  monoxide,  hydrogen  and  oxygen  are  fed  into  the 
combustion  chamber  and  ignited.   The  amount  of  each  reactant 
is  determined  by  the  following  reaction; 


90N2  +  9CO  +  H2  +  502  ■*  9C02  +  H20  +  90N2 


After  the  reactant  products  reach  equilibrium,  nitrogen  is 
mixed  with  the  carbon  dioxide  and  water  in  the  mixing  chamber 
This  yields  a  molar  mixture  of  90%  nitrogen,  1%  water  and 
9%  carbon  dioxide.   This  is  the  composition  required  for 
operation  of  a  gas  dynamic  laser  amplifier.   It  is  assumed 
that  the  ratio  C  /C   of  the  mixture  is  1.4  while  the 
vibrational  mode  of  carbon  dioxide  is  frozen  out  at  T^ . 

The  nozzle  array  is  designed  to  attain  Mach  1.0  at  the 
throat  (region  3  in  Figure  1)  and  Mach  4.0  or  greater  inside 
the  cavity  (region  4  in  Figure  1) .   For  purposes  of  this 
thesis,  all  calculations  have  been  done  assuming  that  a 
Mach  number  of  4.0  exists  inside  the  amplifier  cavity.   The 
technique  used  to  determine  the  area  ratios  is  described  in 
detail  in  reference  [6],   Since  these  calculations  are  well 
documented,  only  some  of  the  final  results  will  be  mentioned. 
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From  the  equilibrium  combustion  temperature,  the  throat  and 
cavity  temperatures  can  be  determined  and  the  small  signal 
gain  can  be  calculated. 

The  small  signal  gain  is  applicable  only  inside  the 
cavity  (region  4  in  Figure  1) .   Assuming  that  the  oscillator 
and  cavity  are  matched  at  a  wavelength  of  10.6  micrometers, 
the  output  beam  intensity  is  given  by  Equation  1. 


I  =  I1  eal  (1) 


a  is  the  gain  coefficient  and  L   is  the  overall  beam  pathlength 
inside  the  cavity.  I   is  found  from  Equation  2. 

I   =  NL  (2) 

N  is  the  number  of  times  that  the  beam  traverses  the  width 
L  of  the  cavity  (see  Figure  1) . 

A  measure  of  the  power  that  can  be  extracted,  defined 
in  Equation  5.23  of  reference  [2]  as  the  limiting  power 
extraction,  is  a  I   where  I   is  defined  by  Equation  3. 

'       o  o  O  J    1 

I   =  k(e3-6l)QVIB  (3) 

(^VT4»£3Eo 

The  maximum  available  energy  in  the  excited  state  of  the 
gases  is  described  in  reference  [2]  and  can  be  calculated 
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with  Equation  4.   Table  II  shows  the  result  of  calculations 
at  incremental  temperatures. 

Emax  e  =  I  3380   1  (4) 

m  e— = 1 

l3 

m  is  defined  as  the   average  mass   of  the   gas  molecules,      m 
is   calculated  with  Equation   5. 

_       *co/co2  +  VV  +  \\ 

m  = (5) 

A 
v 

tj>  is  the  mole  fraction  of  the  subscripted  quantity. 

As  a  basis  for  determining  the  change  in  limiting  power 
output,  small  signal  gain  and  maximum  available  energy,  an 
assumption  was  made  that  all  gases  were  at  a  room  temperature 
of  300  K,  prior  to  combustion.   The  nitrogen  temperature  was 
then  incremented  and  new  calculations  made  at  each  value  of 
nitrogen  temperature.   These  new  calculations  are  compared 
to  those  at  300  K  in  tabular  and  graphical  form. 

B.   SMALL  SIGNAL  GAIN 

Reference  [2]  outlines  procedures  for  determining  the 

small  signal  gain  coefficient.   The  overall  gain  I/I-i  may 

be  determined  from  Equation  1  after  a  has  been  found  from 

Equation  6.   Pages  5-12  of  reference  [2]  shows  the  relation 

between  a  and  a  . 

o 
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-j(j+D  e 
2  r 

n      _  1  Xo   (2J+1)  o    T4    rn(001)-n(100)1  ,    M 

ao  "  *?  ^r  v^7  [~^7 ]  *co2 

(6) 


The  population  inversion  is  determined  by  the  quantity  in 
brackets  in  Equation  6.   Population  inversion  is  related  to 
cavity  and  throat  temperature  by  Equation  7.   It  is  assumed 
that  the  vibrational  temperatures  of  nitrogen  and  carbon 
dioxide  are  frozen  at  T3 .   Appendix  3  is  a  derivation  for 
v  in  terms  of  temperature  and  number  density. 


-e  /t    -e  /t 

rn(001)-n(100)  ,    re  J   J  -  e  x   q.        ._. 
VC02  UVIB 


Here,  the  temperature  T_  is  frozen  at  the  nozzle  throat  and 
the  system  is  operating  so  that  the  threshold  condition  is 
exceeded.   That  is,  63/61  <  To/T4 •   This  condition  insures 
that  the  small  signal  gain  will  be  positive.   Q  .,  is  related 
to  throat  and  cavity  temperatures  by  Equation  8 . 

"61/T4  -1     "G2/T4  -2     "e3/T3  -1 
Qvib  -  d-e  1     4)  X(l~e     4)  Ul-e     3   3)  L  (8) 


C.   MAXIMUM  AVAILABLE  ENERGY 

The  maximum  available  energy  is  dependent  upon  the 

wavelength  of  the  laser  and  the  vibrational  temperature  of 

the  transition  level  for  stimulated  emission.   E    can  be 

max 
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calculated  using  Equation  4  or  Equation  9.   Equation  9  is 

c 

the  same  as  Equation  4  except  that  y  has  been  substituted 

for  v. 


E     -  i  h£  (e  3   3  -  1)  (9) 

max   —  > 

m  A 


D.   HEAT  EXCHANGER 

The  heat  exchanger  is  the  key  to  success  of  this  system. 
It  takes  energy  from  the  exhaust  gases  at  the  diffuser  and 
transfers  it  to  the  nitrogen  gas.   This  process  allows  the 
nitrogen  to  be  heated  to  an  optimum  temperature.   Equation 
13  describes  the  pressure  drop  across  the  exchanger.   Since 
this  system  is  to  operate  at  mean  sea  level  and  higher,  the 
minimum  pressure  on  the  downstream  side  (region  7  of  Figure 
1)  of  the  exchanger  must  be  greater  than  1  atmosphere. 

Assuming  that  the  nitrogen  enters  the  exchanger  at  3  00  K 
and  exits  the  exchanger  at  a  temperature  of  T   (optimal) , 
the  total  heat  transfered  is  defined  in  reference  [3]  as 
Equation  10. 

%    =  ™N  CpN  [TN(H°T  "  VC0LD>1  (10) 

qN  must  equal  the  heat  given  up  by  the  mixture  flowing  through 
the  heat  exchanger.   This  is  described  by  Equation  11. 


q   =  q   =  m   C    [T  (HOT)  -  T  (COLD)]  (11) 

^N   ^m    m  pm   m         m 
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Therefore,  Equations  10  and  11  can  be  combined  to  produce 
Equation  12  which  describes  the  temperature  drop  across  the 
heat  exchanger. 


T5  -  T?  =  (3)  (_P*L)   [TN(HOT  -  TN(COLD)]  (12) 

m    C 

m    pm 

Here.  Tc  is  defined  as  the  upstream  temperature,  T   (hot), 
->  m 

and  T_  is  the  downstream  temperature,  T   (cold) . 

A  more  detailed  analysis  is  in  Appendix  A.   Equation  13 
is  used  in  reference  [4]  to  describe  Ap. 

2 

G  Vl  9      V9  A   V  9      V9 

Ap  =  -~-L   [(k  +l-az)+2(~-l)+f^-  —  -  (1-a  -K  )— ] 
*  zq       L   c  vn       A   v,  e  vn 

1  C   1  1 

(13) 
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III.   CALCULATIONS 

A.   THE  GAS  DYNAMIC  LASER  CAVITY 

Calculations  for  the  combustion  chamber  and  nozzle 
throat  were  computed  by  using  the  NASA  Lewis  Computer 
Program,  reference  [lj .   The  remaining  calculations  were 
S  done  by  hand  using  the  isentropic  tables  in  reference  [6] . 

For  conditions  in  region  4,  shown  in  Figure  1,  Appendix 
A  of  reference  [6]  was  used.   At  M  =  4.0  the  following 
relations  are  found: 


-   Y 


P4  y-1       2       Y  -2 

JL  =    (1    +   JL±  mz)  =    0.6586    x   10 


P  2 

1 


T.  ,   7    -1 

jJ-  =    (1  +  I^i  M  )    =  0.2381 


A      A4 

h- =  a— - —  =  10-72 

A    Athroat 


Conditions  in  region  4  (Figure  1)  can  be  determined  from 
the  following  relationships: 


P4  -2 

P4  =  (P^)P1  =  (°-6586  x  10   )  p!  <14) 

T. 

T4  =  {rf-)^1   =  (0.2381)  T±  (15) 
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A4  =  'a^T1- )Athroat  "  (10-72>  \hroat      (16) 
throat 


One  model  of  the  starting  process  uses  a  normal  shock 
which  develops  inside  the  cavity  (region  4  of  Figure  1) . 
Therefore,  the  diffuser  must  be  able  to  swallow  this  shock 
before  the  system  will  be  stable  enough  to  operate  effec- 
tively.  From  Appendix  B  of  reference  [6]  at  M.  =  4.0,  the 
Mach  number  of  the  flow  behind  the  normal  shock  is  0.4350. 
Before  the  diffuser  can  swallow  the  shock,  the  flow  must  be 
accelerated  to  Mach  1  at  the  diffuser  throat.   Therefore 
from  Appendix  A  of  reference  [6]  at  M  =  .4350,  the  following 
ratio  is  found: 

A     A4 

—  =  — =  1.5  (17) 

AN   Adiffuser 

Diffuser  conditions  change  after  the  shock  has  been  swallowed 
Therefore,  the  only  noteworthy  condition  is  the  fact  that 
the  ratio  of  diffuser  throat  area  to  nozzle  throat  area  can 
now  be  determined. 


'diffuser        A//AN        10.72 


■  . 

=  7.15  (18) 


\hroat   "  S7*5  '  TT 


By  the  methods  described  in  Section  3.4  of  reference  [6], 
the  Mach  number  at  the  diffuser  throat  is  now  3.55  and  the 
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normal  shock  is  positioned  just  past  the  diffuser  throat 
near  region  6  of  Figure  1. 

B.   SMALL  SIGNAL  GAIN  COEFFICIENT 

The  small  signal  gain  coefficient  must  be  determined 
before  the  output  intensity  can  be  calculated.   Combining 
Equations  6,  7 ,  and  8  yields  Equation  19.   Equation  19  is 

explicit  in  terms  of  known  throat  and  cavity  temperatures. 

5  .  . 

The  constant  1.6  x  10   is  the  result  of  inserting 

X   =  10.6  x  10"6m.  J  =  20,  6   =  2.88  K,  v   =  39.5  x  10~14N  JY? 
o  '  '      r  c  v  4 

and  x  =  5.38  sec. 

-239.4        -3380     -1995 

5e4                    e3-e4 
aQ  -    [1.6x10    ][ — __2][ zrggs- — - z^ny ^M ] 

4       -1  4    -2  "}       -1 

(1-e        q    )    ±(l-e        4)    z(l-e        J    ) 

(19) 


Table  1  gives  the  values  of  nitrogen,  throat  and  cavity 
temperatures  used  in  calculating  a   at  the  specified  nitrogen 
temperatures.   Figure  2  is  a  graph  of  a  versus  temperature 
of  the  nitrogen.   It  is  shown  that  an  operating  nitrogen 
temperature  of  approximately  1300  K  will  maximize  the  small 
signal  gain  at  0.773.   See  Table  II  for  values  of  a  at 
specified  nitrogen  temperatures. 
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C.   MAXIMUM  AVAILABLE  ENERGY 

The  maximum  available  energy  can  be  easily  calculated 
by  using  Equation  4  after  m  has  been  determined  by  using 
Equation  5. 


™  -    (.09)  (44)    +    (.01)  (18)    +    (.90)  (28)  .    Q7  v    in-23      gm 

m  -  ^-55 —   4.0/  x  10  — 2, -- 

c   nn        in23  molecul< 

6.02  x  10 

3380 


T 

•1-.    _  /o  o->    ,«5,  ,    3  ,.-1  Joules 

Em=v    (3.83  x  10  ) (e  -  1)    — 

max  Kgm 


These  results  are  listed  in  Table  II  and  Figure  (3)  is  a 
plot  of  E    versus  nitrogen  temperature. 

D.   LIMITED  POWER  EXTRACTION 

The  saturation  beam  intensity,  as  defined  by  Equation  3, 
is  a  valuable  aid  in  approximating  the  limiting  power  extrac- 
tion.  Equation  3  must  be  used  in  conjunction  with  the 
following  equations: 


-e  /t    e  /t    -1    e    e  /t     -i  -1 

ex  =  e  L      q[{e   x      4-l)    +  2  ^±  (e  *  -  1)   ]      (20) 


-e3/T 

3  =  (1  -  e   J   J)  (21) 

-0  J(J+1) 

2 


Z   =  1  _2_i2Jil)_e      4  (22) 

O    4T  xvc  T4/26 
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Combining  Equations  3,  8,  20,  21  and  22  yields  Equation  23. 


"!i        "!i        "e3 

,   _k(e3-fll'   (l-e'VNl-e'V^l-e13)-1         (23) 

°   (!ir1+,.)  1  Xo  2J+1   .-erJ(J^l)A4 

E3  1  4  37  HT  T4/2Sr  e 


The  product  of  Equation  19  and  Equation  23  is  Equation  24. 


m  m 

a   I  9-0  x3  4 

°°   =  [10*         ]  [-i— i]  [2-_ S ]  (24) 

P  2  4  1 

4  e3(4TlP4+T4P4' 


2 

a  I  /P.  is  in  Mks  units.   x.P,  and  t.P.  can  be  found  in 
o  o   4  14      4  4 

Figures  2.8  and  2.10  of  reference  [2].   Table  III  shows  the 

calculated  results  from  Equation  24  and  Figure  4  is  a  plot 

2 

of  a  I  /P„  versus  nitrogen  temperature. 

O  O   4 

E.   HEAT  EXCHANGER 

Using  the  optimum  temperature  as  shown  in  Figure  4 
the  hot  nitrogen  temperature  is  specified  as  2400  K.   This 
gives  a  temperature  difference  in  the  nitrogen  of  2100  K. 
An  energy  balance  will  produce  an  equation  relating  the 
temperature  difference  from  the  hot  to  the  cold  side  of  the 
heat  exchanger.   From  equation  10  and  C   equal  to  0.39, 
the  rate  of  heat  transfer  is 
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•     •  •  RTU 

<3  =  "k  CPN  (iV  =  "fc  (819)  IEi=sec 


Reference  [7]  contains  values  for  C  .   From  Equation  12, 
the  temperature  difference  is 


C  M 

Tc  -  T_  =  (1803)  r^-  z    1617  K 

pm 


A  more  detailed  analysis  of  the  heat  exchanger  is  in 
Appendix  [A] . 

F.   APPROXIMATE  WEIGHT  AND  VOLUME 

Assuming  one-second  bursts  of  energy  in  the  output,  the 
amount  of  fuel  needed  for  any  pre-determined  number  of 
bursts  can  be  calculated  easiest  at  the  nozzle  throat  where 
gas  properties  are  well  defined.   From  reference  [6] ,  the 
mass  flow  rate  is 

m  =  UAp 

where  U  is  the  mass  velocity,  A  is  the  area  of  the  nozzle 
and  p  is  the  density  of  the  gases. 

The  limiting  power  per  unit  mass  is 


a  I 
m 
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where  m  =  (1  sec)  m  .   If  W  increases  by  some  factor  A  and 

m 

m  remains  constant  then  a  I   must  have  increased  by  A  also. 

o  o 

It  is  desired  to  maintain  W  constant  and  adjust  the  mass 
necessary  to  compensate,  for  this  increase  in  output. 


(1  +  A)  a  I 

W(l  +  A)  =  °-° 

m 


a  I     a  I 

oo     oo^^i-   m 
___  -      ■>  m  _  1^ 

m      m  * 
TT+A) 


This  means  that  the  mass  necessary  to  maintain  a  constant 
output  can  be  decreased  by  a  factor  of  m +A\  •   Calculations 
have  shown  that  A  -   0.42  and  the  fuel  necessary  is  only  70% 
of  the  fuel  used  without  a  regenerator.   Therefore  a  savings 
in  fuel  weight  and  volume  of  3  0%  has  been  found. 
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IV.   RECOMMENDATIONS  AND  CONCLUSIONS 

It  has  been  shown  that  pre-heating  the  nitrogen  gas 
prior  to  mixing  it  with  the  combusted  products  of  hydrogen, 

oxygen,  and  carbon  monoxide  increases  the  gain  of  a  carbon 

-1  -1 

dioxide  gas  dynamic  laser  from  0.35  m  '  to  0.77  m   .   The 

limiting  power  extracted  reaches  a  maximum  at  a  temperature 

of  2400  K.   This  produces  an  optimum  operating  temperature 

and  the  limiting  power  extracted  can  be  increased  by  a  factor 

of  1.42.   This  implies  an  increase  of  42%  in  output  power 

while  maintaining  a  constant  mass  flow  rate.   This  means  that 

a  constant  level  of  output  power  can  be  maintained  while 

reducing  the  amount  of  fuels  needed  by  30%.   This  decreases 

the  amount  of  storage  volume  by  30%. 

It  is  recommended  that  a  further  study  be  conducted  to 

determine  the  effects  that  a  heat  exchanger  would  have  on 

other  fuels  to  produce  a  10.6  micrometer  wavelength  with  a 

possible  extension  into  other  wavelengths. 
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TABLE  I 

TEMPERATURES  AND  PRESSURES  OF  SYSTEM 


h  (W 

Tj,  (K) 

T3  (K) 

T4  (K) 

P1=P2(PSIAJ 

P4(PSIA) 

300 

1076 

928 

256 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

5.927 

1000 

6.586 

600 

1303 

1130 

310    i 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

5.927 

1000 

6.586 

900 

1535 

1338 

365 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

5.927 

1000 

6 .  586 

1200 

1774 

1553 

422 

p;nn 

3,293 

600 

3.951 

700 

4.610 

[    800 

5.268 

900 

5.927 

1000 

6.586 

1250 

1815 

1590 

432 

^nn 

3.293 

600 

3.591 

700 

4.610 

800 

5.268 

900 

5.927 

1000 

6.586 

1300 

1855 

1623 

442 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

5.927 

1000 

6.586 

1350 

1895 

1662 

451 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

5.927 

1000 

6.586 

1400 

1935 

1694 

461 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

5.927 

iooo 

6.586 
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TABLE  I  (Continued) 


TN  (K) 

T2  (K) 

T3  (K) 

T4  (K) 

P]=P2(PSIA) 

P4(PSIA) 

1450 

1975 

1730 

470 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

5^927 

1000 

6.586 

1500 

2014 

1772 

480 

5nn 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

5.927 

1000 

6.586 

1800 

2246 

1991 

535 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

5.927 

1000 

6.586 

2100 

2463 

2204 

586 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

90C 

5.927 

1000 

6.586 

2400 

2662 

2404 

634 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

"""'5.927 

1000 

6.586 

2700 

2845 

2590 

677 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.^68 

900 

5.927 

1000 

6.586 

3000 

3017 

2760 

718 

500 

3.293 

600 

3.951 

700 

4.610 

800 

5.268 

900 

5.927 

1000 

6.586 
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TABLE  II 

COMPUTED  VALUES  FOR  Emav  &  a, 

max    c 


TN(K) 

(KJ/KG)    -    Emax    -    (KJ/lbm) 

a0(nrl) 

300 

10.3                      4.67 

0.36 

600 

20.25                     9.12 

0.57 

900 

33.3                      15.1 

0.71 

1200 

49                       22.2 

0.772 

1250 

51.9                      23.6 

0.773 

1300 

54.5                     24.7 

0.772 

1350 

57.7                      26.2 

0.771 

1400 

60.3                     27.4 

0.771 

1450 

63.3                     28.7 

0.769 

1500 

66.8                     30.3 

0.768 

1800 

85.9                      38.9 

0.728 

2100 

105                       47.6 

0.67 

2400 

124                       56.3 

0.61 

2700 

143                       64.9 

0.546 

3000 

159                       72.1 

0.495 
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TABLE   III 

COMPUTED  VALUES  OF  a.Jn/£     (%(— K)  2 

0  °     q       mJ     N/nT 


TN  (K) 

T3  (K) 

T4  (K) 

%W 

300 

928 

256 

0.543 

600 

1130 

310 

0.65 

900 

1338 

365 

0.746 

1200 

1553 

422 

0.827 

1350 

1662 

451 

0.844 

1500 

1772 

480 

0.8834 

1800 

1991 

535 

0.918 

2100 

2204 

586 

0.952 

2400 

2404 

634 

0.963 

2700 

2590 

677 

0.964 

3000 

2760 

718 

0.960 
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Figure  2    SMALL  SIGNAL  GAIN 
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Figure  3    MAXIMUM  AVAILABLE  ENERGY 
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Figure  4    LIMITING  POWER  OUTPUT 
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Figure  5    MIXING  CHAMBER  EQUILIBRIUM  TEMPERATURE 
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APPENDIX  A 
THEORY  OF  GAS-TO-GAS  HEAT  EXCHANGER 

As  stated  on  page  21  of  reference  [4],  the  terms 
involving  exchanger  input  and  output  flow  are  negligible. 
An  inspection  of  Equation  13,  shows  that  the  assumption 


-    f  |-  ^  +  2(!i  -  1)  »  (Ko+l-02)  -  (l-a2-Ke)  !i 
c   1      1  1 

2 

can  be  used  to  develop  Equation  4'.   K  +l-a   is  associated 

2    v2 
with  the  upstream  flow  while  (1-a  -K  ) —  is  associated  with 

e  v1 

the  downstream  flow. 

2 

G  vi     v9  a   v 

AP  =  _Jl  [2(—  -  1)  +  f  £ -]  (4') 

2g     Kv  Ac  vl 


By  definition: 


vm   —2 


Therefore: 


°2vl  V2  f  A  f  A 

AP  =  ^g1  [^{2lT  +   2)  +  <§£-"  2^ 
3    1    c  c 


Now,  the  heat  transfer  area  is  needed.  The  derivation 
is  based  on  the  counterflow  heat  exchanger,  from  reference 
[7] ,  shown  below. 
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TM(COLD) 


The  temperature  variation  as  a  function  of  distance 
through  the  heat  exchanger  is  shown  below. 


T  (K) 


Now: 


TN(COLD) 


DISTANCE 


4  =  hAiT     S 

^       m   sec 


AT   =  logarithmic  mean  temperature  difference. 


h  =  heat  transfer  coefficient. 


A  =  transfer  area 


Before  A  can  be  determined.  AT   is  needed.   Therefore 

m 


d(AT)  =  dT   +  dT._ 
m     N 
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where : 


dT   =  chanqe  in  hot  gas  temperature 
m 


dT   =  change  in  cold  gas  temperature 


Now: 


dq  =  -C  m  dT   =  C„  m.7  dT.T 
^     m  m   m    N  w   N 


Therefore: 


d(AT)  =    dq,   +  ...3SL.. 


m  m 


Also: 


dq  =  hAT  dA 


Therefore : 


d(AT)  _      1        1   »  ,a 
CN  "N    Cm  mm 


Integration  of  this  equation  leads  to 


T7-TN(COLD)  l  1 

ln[Tt.-TM(HOT)  ]  =  hA(CT-mT  "  (TH^ 
5   N  N   N     mm 
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As  shown  previously 


hA  =   q/AT 


and 


q/VTm 


.  WW     ,  Vn[tn(hot)    -tn(cold)] 


AT  AT 

m  m 


Now: 


T„-TM(COLD)  T^T(HOT)     -    TM(COLD)     -    T_    +    T_ 

ir/JN                     ,  N                              JN                                 o              / 

inlT   -TM(HOT)~~J     '  AT 

5      N  m 


and 


[T     (HOT)     -    T     (COLD)     -    T      +    T    ] 
AT        = i__ 

m  T_    -   TrT(COLD) 


In    [= = ] 

5  N(HOT) 


NOW: 


A   = 


m  C    (T,-TJ  pc    Uc   A   C    (Tc    -   T-) 

m  m      5      7'  K5      5      c  m      5  7 


hAT  hAT 

m  m 


And  the  ratio  A/A  can  be  found  as 

c 


[P5U5Cm(T5-T7)1/[hATm] 
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Assuming  a  10%  pressure  drop  across  the  exchanger,  P,- 
equal  to  1.1  atmospheres  and  TV  =  4500 °R,  the  Reynolds 
number  can  be  calculated. 


R  =   S   5   5 
e 


Assuming  a  tube  diameter  of  .1  ft,  Ur  •-  988  (ft/sec)  , 

_  AO    ..-4  SLUG    ,      -  no    1n-7  SLUG 

p  =  3.02  x  10    T  and  u  =  1.09  x  10    kt —  ; 

ft 


R  =  2.73  x  105 
e 


From  Figure  6-9  of  reference  [3] 


hD  9 

-j-£  =  5.75  x  1(T 


k  =  k   .  i/T/T   7 
ref  v    ref 


where  k   ,  and  T   r  are  from  Table  a-5  of  reference  [3] 
ref      ref 


.     /n    nrtrtn.     /4500'    nt;   BTU 
k  =  (0.0337)  ^  =  .05 


hr  ft  °R 


h  =  336  BTU/(hr  ft2  °R) 


p   U   C 

fL.  =  b       b      m  (v   -  T  )  ~  12 
A       "  ._     U5     1}     -    L 
c     hAT 
m 
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v2 
With  these  quantities  and  ^—  equal  to  1.36  and  f  =  0.02 

1  2  4 

^    (from  Figure  6-3  of  reference  [4]  and  G  v  =  1.365  x  10 

lb2/(lbm  ft  sec2) ; 

2 

G  v 

AP  z   ,  •*•  =  0.1  atmospheres  =  212  =£*■ 
^y  ft^ 
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APPENDIX    B 
DERIVATION    OF    v 


From  Equation   3.7   of   reference    [2],    v      is 


l  N.a.C.    =  Nl^a^C^   +   ^2%2CK2   +   *H  0°H  0<TI2ol 


where 


ip  =  the  mole  fraction 


a   is  the  optical  broadening  collisional  cross 
section. 

C.   is  the  relative  mean  speed. 


*C02  =  -09      \   =  -90      *H20  "  -01 


From  pages  3-5  of  reference  [2] 


-14    2 
o    =  1.3  x  10    cm 

-14    2 
a„   =   .87  x  10    cm 

N2 

a„  ^  =   .38  x  10~   cm 
H2° 
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From  page  3.4  of  reference  [2] 


where 


Now: 


c.  =  p~^£ 


y  =  M    for  C02  collisions 

Mco2  Mx 

y  =  for  H20  and  N„  collisions. 

MC02+Mx 


/8   K  T71 
uco  y  it  M  Z-LV14     sic 


/8    R  T„'  M 

F  -     |  i=    44.2\/"^7    -JL- 

N2  K    ■  V      4      sec 

2 
/8   R  t; 

2  I     ir.y 


Therefore 


v      =  N  /T^1  x    10"14[(21)  (.09)  (1.3)    +    (44.2)  (.9)  (.87) 


+    (40)  (.01)  (.38)] 


=   3.95   x   10"13   N  fF^ 
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